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Abstract 

Purely rotational relative equilibria of an ellipsoidal underwater vehi- 
cle occur at nongeneric momentum where the symplectic reduced spaces 
change dimension. The stability these relative equilibria under momen- 
tum changing perturbations is not accessible by Lyapunov functions ob- 
tained from energy and momentum. A blow-up construction transforms 
the stability problem to the analysis symmetry-breaking perturbations of 
Hamiltonian relative equilibria. As such, the stability follows by KAM 
theory rather than energy-momentum confinement. 

Introduction 
The phase space TSE(3) with Lagrangian 

L(A,a,tl,v) = -fi*IQ+ -u*Mw (1) 
approximately models the motion of a neutrally buoyant vehicle submerged in 



an inviscid irrotational fluid (see Leonard [ 1997 1 and the references therein), in 



the case of coincident centers of mass and buoyancy. Here tangent vectors of 
SE(3) are represented by left translation and elements of SE(3) parameterize 
the configurations of the vehicle by embedding a reference vehicle into the fluid. 
I and M are constant, positive definite, 3x3 matrices that can be calculated 
from the shape and mass distribution of the vehicle. This system admits the 
SE(3) symmetry of the left action of SE(3) on itself. 

For an ellipsoidal vehicle with principle axes of inertia along the axes of 
symmetry of the ellipsoid, I = diag(ii, I 2 , h) and M = diag(Mi, M 2 , M 3 ) (i.e. 
I and M are diagonal). If Mi — M2 and I\ = I2 (or similarly if Mi = M3 
and I\ = I3 etc.) then there is a further material symmetry of the system: 
50(2) = { exp(k A (9)} acts as a subgroup of SE(3) by inverse multiplication on 
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the right. In the case that I and M are both constants of the identity then the 
material symmetry is 5*0(3). 

Lie-Poisson reduction yields the Poisson phase space se(3)* = {(7r,p)} where 
7r = IS7 and p — Md. The equations of motion are 

dir „ dp _ 

- = nxQ+pxv, ^= P xf7, (2) 

and by direct substitution, for each a e € M, 

p" s : 7r = a e k, p = 0, 

is an equilibrium of the Poisson reduced systems and hence a relative equilibrium 
of the original system. The generator is 

^3 

so the relative equilibrium corresponds to a stationary vehicle rotating about an 
axis of symmetry which is aligned with the vertical. This article is concerned 
with the stability of these relative equilibria in the case where I3 is not an 
intermediate axis of symmetry, i.e. assuming 1% < I2 < I3 or J3 < I2 < I\. 

The symplectic leaves of se(3)* are as follows. On the complement of p = 
lie the generic symplectic leaves, all diffcomorphic to TS 2 , and which are the 
level sets of the two Casimirs \p\ and ir-p. Nongeneric leaves occur within the set 
p — and are the level sets of the subcasimir |-7r|. Thus the relative equilibria p"' 
correspond to Lyapunov stable equilibria on the (nongeneric) symplectic leaves 
of se(3)* since the energy has a definite critical point when restricted to those 
leaves. Were the symmetry group to be compact this leafwise stability would 
imply stability of the equilibrium modulo the isotropy group of the momentum. 
SE(3), of course, is not compact. The question is whether or not p" e are stable 
under perturbations from nongeneric leaves into nearby generic leaves. 



Leonard and Marsden [ 1997 1 have identified this que s tion as particularly 



delicate, and the theory of Patrick, Roberts, and Wulfi [2002], the sharpest 



possible for the problem of establishing the stability of relative equilibria by 
energy-momentum confinement in the case of noncompact symmetry, corrobo- 
rates that opinion. Patrick et al. separate generators of relative equilibria into 
two complementary classes, tame and wild. The generator of p" " is tame if and 
only if a e — 0, corresponding to a stationary, nonrotating vehicle, in which case 
S!E(3)-stability follows directly, since the energy has zero derivative at p™" and 
has positive definite Hessian there. However, if a e ^ 0, the generator is wild 
and the theory does not imply stability. 

So it is an open question whether the relative equilibria p™ e , a e ^ 0, are 
S'£ ; (3)-stable or not, and the problem appears inaccessible by energy-momentum 
confinement. This is due to the presence of a noncompact symmetry group and 
wild generators. 
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§1 The blow-up construction 



In its essence the stability issue is one of perturbation from a nongeneric sym- 
plectic leaf to nearby, higher dimensional, generic leaves. In order to bridge to 
Hamiltonian perturbation theory, which is usually cast in a setting of a fixed 
canonical phase space, it is natural to begin by normalizing the generic leaves. 
The leaf corresponding to 

\p\ = a, n-p = b, 

for a > is diffcomorphic to 

TS 2 = {{w,w) G M 3 x R 3 : \w\ = l,w w = 0} 

by the map 

p tt ■ p 

W = ]p-\> W = *-^ 

the inverse map being given by, for fixed a > and b G M, 

• b 

p = aw, tt = w H — w. (3) 

Having normalized the symplectic leaves to the constant manifold TS 2 , one seeks 
to extend this to the nongeneric leaves within p = 0, which means extending it 
to a = 0, since p = aw. As it stands (^) is poorly defined for a = 0, but for 
fixed ratios of b/a it is well defined even for arbitrarily small a, suggesting that 
the proper way to approach the nongeneric leaves from generic ones is through 
constant tt ■ p/\p\. Setting 7 = b/a and using the parameters a and 7 instead of 
a and b codes the generic leaves so they fit smoothly into the nongeneric ones, 
thus allowing the possibility of an effective perturbation approach. The map 
p = aw, tt — w + "fw for a = is many-to-one and so the three dimensional set 
of nongeneric leaves p — is "blown-up" by this map to the five dimensional 
set of TS 2 x E = { ((w,w), 7) } . Thus one is led to define the blown-up space of 
se(3)* as 

P = TS 2 x M> x R= {(iu,iu,a,7) : \w\ = l,w- w = 0, a > 0} 
with blow-down map 

p = aw, tt = w + jw 
and corresponding blow-up map, defined on the generic (p 7^ 0) leaves only, 

p . tt-p b 

W = T1 , w = tt--—-p, a=\p\, 7=-. 
\p\ \p\ z a 

The blow-up map is a diffcomorphism from the (open) set of generic leaves to 
the (open) set a > in the blown-up space (the generic sector) , such that each 
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generic leaf is sent to the constant manifold TS 2 . The evolution of the generic 
leaves is transformed to a evolution on TS 2 parameterized by the Casimir values 
a and 7. The blow-down map takes the set a = in the blown-up space 
(the nongeneric sector) to the set of nongeneric leaves, and is many-to-one on 
that sector. Increasing the parameter a from zero corresponds to leaving the 
nongeneric leaves and moving to the generic ones, while 7 parameterizes the 
possible avenues of departure. 

The utility of the blow-up to support perturbation arguments depends on 
whether or not the dynamics of the generic sector can be continued smoothly 
to the nongeneric sector. On the generic sector the vector field that generates 
the dynamics is 

1 dp 1 __i __i , . . 

— — -p X I i = wxl (w + 7UM, 

a at a 

din dw 
~dl ~ 7 ~dt 

t-1 .,-1 dw (4) 

7T X I -K+pXNl p—"/—T 

at 

(w + 710) x I _1 (u) + 710) + aw x M _1 a«j — "fw x I _1 (tt) + jw) 
w x I _1 (ti> + "fw) + a 2 w x M _1 w. 

Obviously this is smooth in a for all a > 0, as required. Dynamics on a = that 
is robust enough to continue through perturbation to small positive a will have 
implications for the original system. By continuity in a, the blown-up vector 
field is a lift by the smooth blow-down map of the vector field for the original 
system, even through the nongeneric sector. Thus the flow on the nongeneric 
sector corresponds through the blow-down map to the flow of the original system 
on the union of the nongeneric leaves. 

Actually, the blow-up has a very transparent reformulation, since P is dif- 
feomorphic to S 2 x R 3 x M> = { (w, n, a) } by the map tt = jw + w. Through 
this diffcomorphism the blow-down map is simply p — aw, which is to say that 
w by itself is enough to desingularize the foliation by symplectic leaves, but not 
enough to normalize the leaves. The blow up map is a proper map since the 
map (a, w) t— > aw is proper. 

Some exploration of the nongeneric sector of the blown-up space may be 
helpful for visualization purposes. For fixed tt the equation 7W + w — tt has 
solution w E S 2 and 7 = ttq ■ w. Consequently the blow-up of the point p = 0, 
7r = tt is in all cases a two sphere. This sphere intersects fixed 7 such that 
I7I < |-7To I in a circle, 7 = ±|7r | in a point, and |-y | > |-7r | not at all. Thus 
departure from the point p — 0, 7r = ttq along 7 > 1 7To | is impossible. The 
nongeneric symplectic symplectic leaf \ir\ — r > blows up to -f 2 + \w\ 2 = r 2 
which is diffeomorphic to S 2 x S 2 , and which for fixed |-y | < r is a circle bundle 
and for I7I = r is a sphere. 



dw 

~dt 

dw 

~dt 
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The Hamiltonian pulls back through the smooth blow-down map to 

H = -t^I" 1 ^ + -^p t MT l p = H° + a 2 H\ 

H° == -(to + 7w) t I~ 1 (to + 7w), H 1 = iw'M -1 ^, 

H is written this way in anticipation of perturbation arguments from a = 
to small nonzero a. The symplectic form Co on the nongeneric sector can be 
calculated from the formula for the coadjoint orbit symplectic forms of SE(3) 
in 



Marsdcn and Ratiu [1994], with the result that 



Co(w, w)((8wi,Swi), (5w 2 ,Sw 2 )) 

= —w ■ (Swi x 5w2 — Sw2 x Swi) — 710 • (5wi x 61112). 

By continuity, the relation = rf-ff persists from a = to a > 0, so the 

vector field (Eh is Hamiltonian at a = with symplectic form Co and Hamiltonian 
_ff°. Thus the evolution on the nongeneric sector is Hamiltonian, in a way that 
smoothly continues the Hamiltonian structure of the generic sector. 

The dynamics on the invariant submanifold p = in the original space P 
admits the subcasimir |tt|. This conserved quantity (conserved on p = only) 
pulls back to a conserved quantity \w + ~/w\ for the nongeneric sector of the 
blow-up space P. Since w ■ w = and \w\ = 1, this gives the conserved quantity 
|w| 2 and hence the conserved quantity /(|to|), where / is any function. The 
Hamiltonian vector field of /(|to| 2 ) is 

dw /'(H) . dw f'(\w\) 

— = - ;[. u wxw, — = -7 - 1 . , u w x w. 
dt \w\ dt \w\ 

Note that rh = w + 71/7 is conserved by these equations, so that 

dw f'(\w\) _ dw f(\w\) _ 

— = rri~ m xw ' 17 = -7 1 ■ 1 m x w, 

dt \w\ dt \w\ 

the solution of which is rotations about m. To normalize the period at 2ir and 
the righthand sense about rh, choose 

to \w\ 



which gives /(|w|) = ~\/^l 2 + \w\ 2 . Thus the nongeneric sector has an addi- 
tional SO (2) symmetry, which acts by 



9 ■ (w, w) = (ex-p(m^O)w, exp(m A 8)w) , m 
and has momentum 

J=-Vl 2 + \w\ 2 . 



2 j. 1^12' 
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This extra 5*0(2) symmetry arises from a subcasimir of the original system. The 
action and the corresponding momentum are defined on the nongeneric sector 
where 7 and w are not both zero. The set where a = 7 = and w — exactly 
corresponds through the blow-down/up to the set where p = and 7r = 0, so 
that the action and its momentum are defined only on an open subset does not 
affect the analysis of the relative equilibria p™ " . 

Here are some aspects of the 50(2) action and its relation to the blow- 
down/up map. 

1. The action is free except on the set w = 0, which is a 2-sphere of fixed 
points. This two sphere is also the level — |— y | of the momentum J, is 
a symplectic submanifold of P, and as such is equal to its own singular 
reduction. 

2. The orbit relation of the action together with the parameter 7 exactly 
absorb the additional phase space from blowing up the nongeneric leaves. 
Indeed, for fixed 7 the blow-down map is a quotient map for the action, 
and the orbit space is therefore smooth, irrespective of the fact that the 
action is not free. 

3. The blow-down map restricts to a quotient map for the (singular or non- 
singular) symplectic reduced space associated to the J = fi level set. As 
such this reduced space is symplectomorphic to the nongeneric leaf p = 0, 

M = -A- 

Only the verification of the third item in the nonsingular (fi < —\"f\) case is 
troublesome. For that, it is easily verified that the map ir = -fw + w is a quotient 
map for the 50(2) action on J _1 (/t) which has image the sphere TS^^ = 
{-7T : \ir\ = —fi}- To pull down the symplectic form Co by that map, first let 
(it, 5-Ki) G TS 2 ^^, i = 1,2, and seek (w, w, Swi, Svbi) such that 

\w\ = 1, w ■ w = 0, w ■ 5wi = 0, 5u>i ■ w + w ■ Swi = 0, 



~~Vl 2 + IH 2 = Aj fiWi ■ W = 0, 7T = JW + W, STTi = jSlVi + Swi. 

To solve these equations choose and w such that w ■ ir = 7 and set w = n — jw. 
Expanding Su>i and Swi in the basis w, w, w x w gives 

it) ' Stt ( uj x 7T ) ■ Sir ' 

5wi = -— l -(n~jw), 5wi = (w ■ 5-Ki)w H — ^ioxit. 

[i z — Y [i A — j z 



G 



Substitution into Co then gives 



u>(w, w)((Swi,8wi), (Sw2,Sw 2 )) 

2 (^(Stti ■ w) ((tt X 5tt2) ■ w) — (8tt2 ■ tf) ((tt X fori) • 



fi 2 -T 



(w X (w X 7f)) • (5tTi X (57T 2 ) 



/t 2 — 7 2 



((w ■ tt)w — 7r) • (fori x for 2 ) 



(w ■ TT)- — |2"7T — 7T ) • (fori X foro,) 



"2 2 



1 / 7 2 \ 

= To o ^2 ~ 1 ^ ' (^1 X 571 " 2 ) 

M -7 \M / 

kr 

This last expression is the symplectic form on the nongeneric leaf | vr j = —fx, as 
required. 

In short, the 50(2) symmetry arises and exactly absorbs the additional 
dimensions resulting from the blow-up construction. The symplectic reductions 
of the nongeneric sector by this symmetry exactly coincide with the original 
system restricted to the nongeneric symplectic leaves of the phase space P. 

The pull-back of the relative equilibria p" " by the blow-down map is the set 
of (w, w, a, 7) such that 

p = aw = 0, tt = a e k = w + jw. 

Since \w\ = 1 the first equation is equivalent to a = (the relative equilibria are 
of course in the nongeneric sector), and dotting the second with w shows it is 
equivalent to 7 = a e k- w and w = a e k — jw. Since 7 2 + \w\ 2 — a e 2 and a e 7^ 0, 
all these solutions are in are within the open set where the 50(2) action is 
defined (i.e. where 7 and w are not both zero). Thus the relative equilibria c 
blow up to 

p"° : \w\ = 1, w = a e k — a e (k ■ w)w, a = 0, 7 = a e k ■ w. (5) 

Since w is unconstrained in (|^), except for the first equation, is diffeomorphic 
to S 2 . Since m = k on p"", the 50(2) action on p"" is by rotation of the 
pair (w,w) about k. By substitution of pi" into (Q), each point of p"" is a 
relative equilibrium for the 50(2) symmetry, except for the the two points 
w = ±k, w = 0, a = 0,7 = ±a e , which are equilibria that reside at singular 
points of the action. Each relative equilibrium in j3" e has the same generator, 
namely 

£- = -f (6) 
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and by substitution into J, the same momentum, namely 

Ae C = -«e- 

Fixing 7, which means fixing a parameter, a e = 7/ cos^, where is the 
angle between k and in. Thus for fixed 7 there are 2-submanifolds of relative 
equilibria, as a e is varied, as expected for an SO (2) symmetric Hamiltonian 
system. Along those submanifolds there is the momentum- generator relation 

t* = jK% (7) 
which will give a crucial component in the KAM twist condition to follow. 



§2 Normal forms in the blown-up system 

Proposition |Tj below relates the stability of the relative equilibrium p" ' to the 
stability of its blow-up p" e . 

Proposition 1 Suppose that, for some fixed a e , p"" is stable for the How F t on 
P, in the sense that for all neighborhoods U ofp" s there is a neighborhood V of 
such that F t (p) £ U for allp £ V. Then is a stable relative equilibrium. 

Proof. Suppose U is a neighborhood of p"" . U pulls back by the blow-down 
map to an open neighborhood U of p™ e . Let V be a neighborhood as in the 
statement of the proposition. Then it suffices to show that V pushes forward 
by the blow-down map to a neighborhood of p™ " . But this follows since the 
blow-down map is proper. ■ 

In particular, if all of the relative equilibria and both equilibria in p" ' are sta- 
ble under perturbation both within the phase space TS 2 and in the parameters 
a and 7, then the original relative equilibrium " is stable. When a is perturbed 
away from this is a SO (2) symmetry breaking perturbation. As dimT5 2 = 4, 
the blown-up system is integrable when a — and hence the stability issue is 
one of the stability of periodic orbits of a near integrable Hamiltonian system. 

Assume, without loss of generality, that a e > 0. Since the 50(2) symmetry 
on " is by rotation about k, it suffices to consider the stability of orbits in 
p" B emanating from points p""' e obtained by substituting w — sin^i + cos 0k 
into (|), for 6 € [0,tt]. 



§2.1 Normal form for the relative equilibria 

Consider first the relative equilibria in p""'®; i.e., exclude the equilibria corre- 
sponding to 9 = and 9 = n. The argument proceeds by adapting and incre- 
mentally refining, to the order required for the stability analysis, the normal 
form near relative equilibria developed in Patrick 1 1995| . 

Below 0(x;y) k , x £ R™, y £ R m will denote the set of smooth y de- 
pendent functions such that 0(x;y)/\x\ k is bounded near 0. The product 
0(x;y) h O(x';y') k denotes the set of finite sums of products of elements of 
0(x;y) k and 0(x'; y' f . 
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§2.1.1 Initial normal form. This is constructed from the linearization of the 
relative equilibrium, which means the linearization at p""' e of Hamiltonian vec- 



tor field Xfi a where 



The characteristic polynomial of the linearization is x ^> x 2 (x 2 + uj e 2 ), where 



-W(i-*)(5-5> 

For later convenience define uj e positive if -Z3 > I\ and 73 > I2 and negative if 
I3 < I\ and I 3 < I 2 . The linearization has a and ±iuj e generalized eigenspaces, 
both of dimension 2. Introducing the parameter 

h(h-h) 

~ hih-hV 

the vectors 
V! = —— \ cos6» aesin 2 ^ Ol 

V2 = — [ cos 9 — sin 9 a e sin 2 6 a e sin cos 9 1 



u 3 = sin6» [0100 -a e cos6> ] 

v 4 = — -[ — cos 2 9 cos 9 sin 9 a e cos 3 9 — a e sin 0(1 + cos 2 9) 1 

a e smfl 



form a basis of T~ ac ,eS 2 which satisfies the following: 



1. the basis is symplectically canonical, so that the symplectic form with 
respect to it is 



10 
-10 
1 
0-10 



2. with respect to the basis the derivative of the momentum is 

dj(p*' e ) = [ 1 ] ; 

3. the third basis vector w 3 is the infinitesimal generator action corresponding 
to 1 e so(2); 
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4. the first two vectors v\ , V2 span the uj e generalized eigenspace and the last 
two V3, W4 span the generalized eigenspace of the linearization; 



5. the linearization of the relative equilibrium is 



dx 60 (p^n = 



where 


































K e 















(9) 



(10) 



Consequently, the basis effects a Witt-Moncrief decomposition 
T.^,eP = Nx so(2) ® so(2)* = span(ui,v 2 ) ®Rv 3 Q 



iVi. 



Here N\ , the symplectic normal, may be identified with the tangent space to the 
symplectic reduced space at p™" f° r the 50(2) action. The appearance of the 
nilpotent part of the linearization is the foundational element of Patrick [1995]. 
The value of n e coincides with the derivative d£^^ /djl^" from (f7), as predicted 
by the general theory. 

The initial normal form can now transcribed from the data above, and is 



H 



1 



(q 2 +P 2 ) t,^ 2 +R{q,P, v) + —H\q,p,<p,v), 



(11) 



R = 0{q,p, vf 



on the product of K 2 x T* 5*0(2) = { (q,p), ((p, v) } with the product symplectic 
form dq A dp + dip A dv , with SO (2) acting by lifts of its left action on itself, and 
with the momentum map v — a e . The transcription is that there is an 50(2) in- 
tertwining symplectic diffeomorphism from a neighborhood of the 50(2) orbit 
of p" c to a neighborhood of times the zero section of T*S0(2) which 

1. sends the relative equilibrium p^" 6 top = q = i' = (p = 0. 

2. has derivative at p" ' the identity map with respect to the basis and the 
standard basis of R 2 x T*50(2) ; 

3. intertwines the momentum maps J and v — a e \ 

Thus the transription is strucure preserving in that it is symplectic and it pre- 
serves the 50(2) symmetry and momentum, so the blown- up system near the 
group orbit of the relative equilibrium p""' 6 can be replaced by the entirely 
equivalent system (|IT|) near q = p = v = 0. 
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$2.1.2 Elimination of qO(v) 2 , pO(v) 2 , and (q 2 — p 2 )v. The remainder term 
of ( pi] ) can be expanded as 

R =cx{v)q + c 2 (v)p + c 3 v(q 2 - p 2 ) + c A v(q 2 + p 2 ) 

+ 0{q,pf + 0[yf + 0(q,p) 2 0(q,p, v) 2 , 

where c\(y) = 0(v) 2 ,ci{v) = 0(v) 2 , and C3, C4 are constants. The transforma- 
tion 

c% pv dci 

q = q-\ , ip = (p-\ — , 

LU e u e dv 

suggested by completing the square in \io e q 2 + ciqv 2 , is structure preserving 
and changes the Hamiltonian to the same form but without terms of the form 
qO(v) 2 . Similarly one eliminates pO(y) 2 . The transformation 

*=wr *= f ^ /h ^(tt£) 4 ' (12) 

takes the fragment ^f{q 2 + p 2 ) + c 3 {q 2 — p 2 )v to 

( Y + C3 ) q2 + ( Y - C3 ) p2 = Y (f+f) + eo{v)2 +p 2 °^ 2 > 

while the symplectic form becomes 

/' 

dq A dp + df A dv =dq A dp + dip A dv + — (qdp + pdq) A dv 

( f \ 
=dq A dp + d I ip + ~j1V J A dv. 



Adjoining (p = tp + {f'/f)qp to ( |12[ ) gives a structure preserving symplectic 
transformation that eliminates the term C3(q 2 — p 2 ). Thus, without loss of 
generality, 

R = c A {q 2 +p 2 ) + 0(q lP f + 0{vf + 0(q,p) 2 0(q,p, v) 2 . (13) 

$2.1.3 Normal form for the rigid body. We will require the first two terms of 
the normal form corresponding to the equilibrium it = a e k of the blown-up 
system reduced by its 50(2) symmetry, i.e. the symplectic reduced spaces of 
the rigid body |-7rl 7T. The map 

tt = ((a e - ±(Q 2 + P 2 )) *P, (a e - \{Q 2 + P 2 ))^Q, a e -\{Q 2 + P 2 ) 

is a symplectic chart on the reduced space |7r| = a e and in these coordinates, 
the Hamiltonian becomes, up to a constant, 
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Action-angle variables for the linearized flow are 

Q = V2IDi sinijj, P = V2ID~i cos ip, 
and the Hamiltonian is then 

-■k^tt = ujJ - — (D? sin 2 ijj + D-? cos 2 ip)I 2 . 
2 2a e 

By averaging over tp, 

^Itt = u; e I+^v e I 2 + 0(Q,P) 5 (14) 

where 

§2.1.4 Matching and normalizing the reduced spaces at p"" e . For a = the 
symplectic reduced space through p™ e ' 8 of the blown-up system is the \ir\ = a e 
symplectic reduced space of the rigid body ^^Iit. For a — the symplectic 
reduced space of the normal form ( p"l| ) through q=p = tp = v = 0isM. 2 
with symplectic form dq A dp and Hamiltonian H\ u =$. Since the intertwining 
map between the blown-up system and the normal form is structure preserv- 
ing, it descends to symplectomorphisms of reductions of the these two systems. 
Consequently, by symplectomorphism on {q, p) only, the normal form Hamilto- 
nian ([ll]) atg = p = ^ = can be equated to the rigid body normal form (|l4|) , 
after which the normal is correct to fourth order in pure q and p and 

R = 0(q,p) 5 + is0(q,p) 2 0(q,p,v). 

§2.1.5 Refinement by matching the normal forms and generators along the rel- 
ative equilibria near p""' e . Advantage may be obtained by comparing reduced 
normal along the relative equilibria p^+ z > e as z varies. These relative equilibria 
occur (for both systems) at momentum — (a e + z). For the rigid body the only 
z dependent adjustment is in the a e dependence of the linearized frequency, 
which becomes io e (a e + z)/a e , so the normal form is 

(w. + ^j I+^v e I 2 + 0(q,p;zf. (16) 

For ( |ll"|) it is the normal form of the reduction at v = — z, so it is the normal 
form of the Hamiltonian 

uj e I + ^v e I 2 - 2c A Iz + 0(q,p) 2 0(q,p, z) 2 

which is 

(uj e -2c4Z + 0(z) 2 )l + 0(q,p;z) 3 . (17) 
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Comparison of (|T^) and (|T^) at first order in / gives a crucial fact: 

We 

2c7 e ' 



U! e 

c 4 



Also, the SO(2) generator of the blown-up system at p"' +z ' e , which is ~(a e + 
z)/Is, and the 5*0(2) generator of system (ll]) at the relative equilibrium q = 
p — are the same. Equating these gives 



+n e u+ — 



= - — (a e + z) = C* - K e z, 

q=p=0 i3 



v— — z 



which means that R has no pure v terms. Particularly, the 0(v) 3 term in ( p"3| ) 
is zero. 

$2.1.6 Symmetry breaking term. The transcription to the initial normal form 
is known to first order since it has derivative the identity map along the 50(2) 
orbit of p^"' 6 . Consequently, H 1 can be calculated to first order by substitution 
of 

w = exp(ipk)P w [p"'" e + (qvi +pv 2 + Wi)) 
into w t 'M.~ 1 w, where P w {w, w) = w. 

§2.1.7 Altogether. Putting all the foregoing together, the normal form is 



H =loJ + \v e I 2 + frv + \n e v 2 - *±Iv 
2 I a e 

+ 0{q,pf + vO(q,p) 2 0(q,p,v) 

+ a 2 H lfi (q,p, ip, v) + a 2 H 1A (q,p, <p, v) + a 2 0{q,p, v\ ip) 2 , 



(18) 



where 



and 



H 1 ' =Ml^l S in 2 9 cos 2 p (19) 
2M X M 2 



frll _ (M 2 - M x ) sin 26 



D^qsvn2p + D 4 pcos2</? 



AMxM 2 ^/a, 

cos 2 9/1 2 1 2 1 
— — cos tp + — — sin ip — — — 

Mi Y M 2 M 3 



sin 20 


( 2 


1 






\M 3 


Mi 


M 2 



The details of the symmetry breaking term H 1 are not required for the stability 
analysis and are displayed here for the sake of completeness. The functional 
form of H 1 ' 1 depends on the choice of the basis Vi and further normalization or 
analysis would be required to extract information from it. 
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§2.2 Normal form for the equilibria 

There remains to consider the two equilibria p" e, ° and p"' ,7V corresponding to 
w = k and w = — k, respectively. These equilibria are fixed points of the action 
of 50(2) and the analysis requires a transparent extension of the normal form 
in Patrick| ]1995 | to equilibria which have 50(2) isotropy. 



It suffices to consider p™*" ; the case of p""' 71 ' is similar. There is a one 
parameter family of possible linearizations of the equilibrium, namely the lin- 
earizations at p" c, ° of the Hamiltonian vector fields X^ where — A J. These 
linearizations have characteristic polynomials 

x i > (x 2 + (a e + XI 3 ) 2 )(x 2 +uj e 2 ). 

Choosing A = —a e /I 3 gives the largest possible null space and therefore the 
largest number of intrinsically defined higher order terms. The vectors 

i>i,o = [ 1 ] , v 2 ,o = ^-^= [ 1 ] 



i>3 0=— !=[0 1 -a e 1 w 4 o=^=[-l a e 1 , 



form a basis of X^a e .o5 2 which is symplectically canonical and such that the 
linearization dX H o{p"" , °) is the linearization @ with n e replaced by zero. The 
first two vectors span the tangent space to the (singular) reduced space through 
p" e '°, which is the | -zr j = a e symplectic reduced space of the rigid body ^tt*!^ 1 ^. 

Higher order terms of the Taylor expansion of X H o are intrinsically polyno- 
mials on kerdXj^o (pg e '°). Alternately one can compute the higher order terms 
of the Taylor expansion of the Hamiltonian on the null space. Letting (x, y) be 
the coordinates on ker dX H o(p" s '°) indicated by the last of the two basis vectors 
above, the Hamiltonian on the null space is easily computed to be 

-L(x 2 + y 2 ) 2 + 0(x, yf = ^(x 2 + y 2 ) 2 + 0(x, y) 5 . 

8-Z3 8 

The initial normal form, obtained by the Equivariant Darboux Theorem, is the 
Hamiltonian 

£j W e / 2 1 2\ a e ( 2 , 2\ , K e I 2 , 2\2 

+ 0(q,p)0{q,p, x, y) 2 + 0(x, yf + a 2 H l {q,p, x, y) 

on the phase space K 2 x I 2 = { (q,p), {x, y)}, with symplectic form dq A dp + 
dx A dy, with 50(2) acting by counterclockwise rotation on (x,y), and with 
the momentum mapping — \{x 2 + y 2 ). The transcription is by local sympl ectic 
diffeomorphism with analogous properties to tho se stated in Section 



Manipulations similar to those in Section 2.1 are required, as follows: 



2.1.1 
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1. Linear terms in q of the form qO(x, y) 2 can be removed as in Section 2.1.2 
and similarly linear terms in p, x, and y. Quartic terms of the form 
0(q,p) 2 0(x,y) 2 must by 50(2) invariance be in (x 2 + y 2 )0(q 7 p), and so 
can be written as sums of (q 2 + p 2 )(x 2 + y 2 ) and (q 2 — p 2 )(x 2 + y 2 ), and 
the latter kind removed, as in Section p. 1.2 . 



2. Pure q and p terms up to order 4 can be found by matching the reduced 
system of the initial normal form to rigid body reduced spaces. 

3. By Items 1 and 2 all terms up to and including order 4 are removed or 
calculated, except for the coefficient of the term (q 2 + p 2 )(x 2 +y 2 )- This 
can be found by matching normal forms along the equilibria q = p = 
(which are fixed points of the action of 50(2)), and the resulting term is 

-itS<i 2 +p 2 )( x2 + y 2 )- 

4. The remainder after all of that, having no terms linear in any variable, 
and being at least degree 5, is of the form 0(q,p,x,y) 2 0(q,p,x,y) 3 , and 
is 50(2) invariant. 



5. The symmetry breaking term H 1 can be calculated as in Section 2.1.6 by 
substituting 

w = P w (p* e '° + qvi t o + pv 2 ,o + xv 3:0 + yv 4 ,o) 
into w t 'M.~ 1 w, and keeping the leading terms, which are order 2. 
Altogether, the normal form is 



H =uo e I+ \vj 2 +^v+ \n e v 2 - —Iv + 0{q,p,x,y) 2 0{q,p 7 x,y) 3 



a 2 



(20) 



t y-ff 1 ' 1 ^,??,^,?/) + a 2 0(q,p,x,y) 3 

where / = \ {q 2 + p 2 ) and v = \ {x 2 + y 2 ), and 

£1,1 _ (Mz-M^D-ip-y) 2 (M 3 -M 2 )(Diq + x) 2 
~ a e MiM 3 a e M 2 M 3 

§3 Stability 

The rescaling 

/ = a 2c I, v = a 2c v 

is symplectic with multiplier a 2c . Substituting into ( |i"8| ) and dropping the tildes 
gives 

-^r-H = lu p I + \vj 2 a 2c + &'v + \-K e u 2 a 2c - loJvcl 2c + 0(a; q,p, (p, vf c 
a zc 2 2 

+ a 2 - 2c H 1 > + O{a;q,p, <p, vf- c 
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Matching the exponents of a in first nontrivial terms of the integrable part, 
i.e. ^K e v 2 a 2c and ^v e I 2 a 2c , with the first term of the nonintegrable part, gives 
2c = 2 — 2c, or c = ~ . After putting e = y/a, and disposing the factor 1 /a 2c of 
77, which merely reparameterizes time, one has 

77 = cuj + £°v+ ( 1 -k c u 2 - ^-Iu + \v e I 2 + H 1 ' ) e 2 + 0(e) 3 , (21) 
\2 a e 2 ) 

where the dependence of O(e) 3 on all of q, p, <p, and v has been notationally 
suppressed. For e = the Hamiltonian ( |2l"| ) has a periodic orbit cylinder by 
varying tp and v with 7 = 0. The orbit I = v = corresponds to the relative 
equilibrium p""' 6 . 

For determining stability it suffices to approximate the Poincare map for the 
orbit corresponding to v = in the zero energy level. Solving ( pl[ ) for v when 
77 = gives 



u e 7 

V = ; 

and the equations of motion for ( ^l| ) are 



0( £ ) 2 



dip , { r W e \ 2 , ^m3 °^ \3 



(22) 

, , -, 2 + 0^ ;i 

e \ a e / rft dip 

On the zero energy level the equation for the evolution of tp is 

ffrfr _ *.(«.« e+e) Je2 + 0(£) 3, (23) 

With an initial condition 7 = 7q, the second equation of ( p2| ) gives 7 = 7 +O(e) 3 , 
so the return time of the Poincare map is, from (^3|), 

T EE -%L ( 1 + i^e+&M e2 \ 3 

Solving the first two equations of (|22] ) over this period, dropping the subscript 
for the initial conditions, and using the relation £" e = — n e a e to eliminate a e 
gives 

7' = 7 + 0(e) 3 

^ = * - - e2 J^ ( e) 2 " ^e 2 )7 + 0(ef (24) 



This is of the form ( Meyer and Hal] [ 1991 1 , Theorem 2, page 231), namely 

(J, V) -»(/', V') by 

7' = 7 + e f -+ s c(7>,e), 

^ - 4> + u + e s h(I) + e s+r d(I, e) 
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Figure 1: Numerical verification of the twist predicted by fl24J). Leftmost: the curved leading 
edge of the numerically computed Poincare map indicates a twist map by visibly showing 
faster rotation as / increases. The twist overlays a constant rotation (in I) which is caused by 
high order terms in a and decreases as a decreases. Center: the rotation angle per iteration of 
the Poincare map on the right as a function of /. The slope corresponds to the twist predicted 
by (E3). Right: reciprocal of the twist for I3 = 1, I2 = I3 — Ii as I\ ranges from to 1 
compared to the parabola predicted by (E4j). 



with r = 1, s = 2, and 



h(I) 



2n 



The twist condition dh/dl ^ is 

which, after substituting (0). (p|), rtToh, and 



Ml?*) 2 -KeW a )J. 



2 
h 



1 

h 



1 

% 
1 



hh 



^0, 

5|), i! 
1 

h 



(25) 



1 

h 



l 

7 2 



/ 3 -x(/i + / 2 ) ^0 



This is certainly true if ^3 is not between I\ and I2 

As for the equilibria p""'° of Section 2.2, Arnold's Stability Theorem (Meyer 



and Hal] [ 1991 1 , Theorem 1, page 235) together with the normal form (20) imply 
that the equilibrium is stable when a — if 



2 a e 



U — — LU e 



K e UJ e 2 ) 7^ 0, 



which is the same as the twist condition (|25|). Stability follows for sufficiently 
small nonzero a since a contributes continuously. 

Thus p" c are stable as equilibria on the Poisson reduced space {(-7r,p)}. 
From Patrick, Roberts, and Wulfl [ 2002 1 , and since the momentum at has 
zero translational part and rotational part parallel to k, this implies 50(2) x C 
stability of the original equilibrium, where C is any cone about k. 

The foregoing sort of analysis will always lead to some stability condition, 
irrespective of possible errors in the derivation, so it necessary to check (E4j) by 
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comparing it with numerically generated Poincare maps. Substitution of e = 1 
after truncation of 0(e) 3 into ( pi| ) gives the leading behavior of the Poincare 
map when a, v, and /, are of comparable order (they are all order e 2 ). The 
Poincare map is determined to leading order in / by the first order twist term 
h(I) in when the zero order term 2wu} e /^"" has a vanishing effect (i.e. is a 
multiple of 2n). This happens to occur when I\ + h = h, as is easily verified. 
After substitution of 1% = I3 — Ii , the first order twist is 



whereupon 

/ 1 



h{I) ttL 



h(h-h), 



which is a parabola in I1/I3. As can be seen in Figure |], this compares well 
with numerical integrations of the original (as opposed to the blown-up) system. 



Summary 



The following theorem has been proved. 

Theorem 1 Within the context of the Lagrangian system the motion of 
an underwater ellipsoid rotating about a long or short principle axis of inertia 
is stable modulo SO(2) x C where 50(2) acts around the rotation axis and C 
is any cone containing that axis. 

This theorem follows from KAM confinement after a blow-up construction and 
normal form analysis rather than confinement by Lyapunov functions derived 
from energy and momentum. 
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